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We extend our previous study of taste-symmetry breaking using HYP-smeared staggered fermions 
in two ways. First, we improve the statistics of a comparison of unimproved and HYP-smeared 
staggered fermions on quenched lattices with a w 0. 1 fm. This allows us to obtain a signal for 
all pion tastes, rather than just a subset, and thus to make a complete comparison. In addition, 
it allows us to differentiate between wall and local sources. Second, we compare HYP-smeared 
valence quarks to asqtad valence quarks on 2+1 flavor unquenched MILC lattices. We find that 
taste breaking is substantially reduced by HYP-smearing, bringing the size of this discretization 
effect (which is the dominant such effect with staggered fermions) down to the size expected 
generically for any fermion type. 
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1. Introduction 

Unimproved staggered fermions suffer from very large discretization errors, particularly those 
that involve the breaking of taste symmetry Reducing these errors is crucial for practical appli- 
cations, and one popular choice is the asqtad action [|l|] which has been extensively used by the 
MILC collaboration. This is a tree-level 0(a 2 ) improved action, and has been found to signif- 
icantly reduce taste breaking. Nevertheless, taste-breaking remains the dominant discretization 
error, and further improvement is desirable. This can be accomplished by the use of multiple levels 
of link-smearing (with reunitarization at each stage) [Q,|3]]. A tree-level 0{a 2 ) improved action with 
multi-level smearing is the HISQ action [^J — and this is the present staggered action of choice for 
charm quarks. For light quarks, however, where a 2 m 2 q effects are likely small, we are pursuing a 
simpler option. This is the "HYP action", in which one uses HYP-smeared links (which have three 
levels of smearing) but otherwise keeps the unimproved staggered action. This action is not 
fully 0{a 2 ) improved, but our hope is that it is effective at reducing the dominant taste-breaking 
discretization error down to the generic size, namely <3 2 Aq CD . 

Last year we presented our first results from a comparative study of taste-breaking with unim- 
proved, HYP-smeared and asqtad valence fermions We used the splitting in the pion multiplet 
as a non-perturbative measure of taste breaking. We found that for quenched configurations with 
a pa 0. 1 fm there was a dramatic reduction in taste-breaking. Indeed, our results indicated that, 
rather than treating the splittings as a leading order effect in a joint chiral-continuum expansion, 

1. e. a 2 ~ p 2 , they should be treated as a next-to-leading-order effect, i.e. a 2 ~ p . This would 
substantially simplify the formulae of staggered chiral perturbation theory [|[ [J. 

During the last year we have improved the statistics of the quenched calculation, allowing a 
more thorough calculation of the pion spectrum, now including all tastes, although the overall con- 
clusions are unchanged. More importantly, we have calculated the pion spectrum with the valence 
HYP action on unquenched MILC "coarse" lattices {a pa 0.125 fm), allowing a direct comparison 
with results using valence asqtad action, and also allowing us to determine the power-counting 
appropriate for the HYP action. 

We do not repeat here the technical discussion concerning the types of sources and sinks that 
we use, refering the interested reader to last year's talk [Qj. We focus instead on the new numerical 
results and the corresponding conclusions. Further details will be given in an upcoming publication. 

2. Update on quenched results 

We use quenched lattices of size 16 3 x 64 generated with the Wilson gauge action at coupling 
/3 = 6, so that I /a pa 1.95 GeV. Compared to last year, we have increased our statistics from 218 
to 370 configurations. For unimproved fermions we use quark masses am q = 0.005, 0.01, 0.015, 
0.02, 0.025 and 0.03, while for HYP fermions we use am q = 0.01, 0.02, 0.03, 0.04, and 0.05. 
Despite appearances, these ranges corresponds to heavier physical masses for unimproved than for 
HYP-smeared fermions. This is because Z m pa 2.5 for the former, while Z m pa 1 for the latter. For 
example, the physical strange quark mass is am^ hys pa 0.025 for unimproved fermions and 0.052 
with HYP smearing. 
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Figure 1: (am n ) 2 vs. am q for unimproved staggered fermions with cubic wall sources on quenched lattices. 
Left panel: LT tastes; Right panel: NLT tastes. Linear chiral extrapolations are shown. 

The increase in statistics allows us to calculate the masses of pions of all tastes. Last year 
we had results only for the "LT" (local in time) tastes, i.e. §5, §4, ^5 and ^4. For pseudoscalar 
spin, 75 = 75, these are the tastes created by an operator that is local in time. This year we have 
signals also for the "NLT" (non-local in time) tastes 1, ^, §45 and £y. These are produced by our 
sources (which reside on a single timeslice) by their coupling to the axial current, 75 = 7475, which 
is suppressed relative to the production of LT tastes. 

As described last year, we use two choices of sink operators, and find them to be indistin- 
guishable. All results presented here are obtained with the "Golterman" [Q] sinks. We also use 
two sources — an extended "cubic wall source" and a local "cubic U(l) source". We find that both 
work comparably well for LT pions, while the cubic wall source leads to smaller errors for the NLT 
pions, although the improvement is smaller for the HYP action than for the unimproved action. 

In Figs. |] and || we show, respectively, the pion spectrum for unimproved and HYP-smeared 
staggered fermions. The dramatic reduction in the breaking of taste symmetry observed last year 
is seen to carry over to the NLT tastes. For all except the taste singlet, this result is expected 
because of the predicted approximate SO(4) symmetry, which holds at leading order in staggered 
chiral perturbation theory []|]. The SO(4) symmetry combines tastes (NLT) and ^4 (LT) into a 
single multiplet with taste and similarly for the SO(4) irreps with tastes £5^ and ^ v . Of the 
NLT tastes, only the taste singlet is unpaired, and thus the fact that its mass is very close to those 
of the other tastes does not follow from SO(4) symmetry, but rather is an additional, independent 
indication of the smallness of taste breaking with HYP fermions. 

In Table gj, we list the mass-squared splittings after linear extrapolation to the chiral limit. 
These are proportional to a 2 in staggered chiral perturbation theory. Note that knowledge of Z m 
is not needed in order to compare these quantities between the two types of fermion. Our main 
conclusion is that taste breaking is reduced by an order of magnitude by HYP smearing. In fact, 
only for the LT tastes (the left columns) are the splittings significantly different from zero with 
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HYP smearing. For the NLT tastes, even the ordering of states is not clear. We also note that there 
is evidence for SO(4) symmetry breaking for unimproved fermions at the 2 — 3a level, with the 
NLT tastes systematically lower than their LT partners. We are investigating this issue further. 



Taste (F) 


Af (unimproved) 


Af (HYP) 


Taste (F) 


Af (unimproved) 


Af (HYP) 




0.0230(09) 


0.0022(09) 




0.0200(34) 


-0.0012(29) 


£» 


0.0360(24) 


0.0044(10) 




0.0245(40) 


0.0003(31) 


^4 


0.0447(27) 


0.0059(22) 


& 


0.0342(42) 


0.0022(30) 








I 


0.0432(44) 


0.0040(30) 



Table 1: Comparison of taste symmetry breaking with unimproved and HYP-smeared fermions on quenched 
lattices with cubic wall sources. Here Af = [am Jl (F)] 2 — [am x (^)] 2 , exttapolated to the chiral limit. 

In conclusion, we confirm and extend the findings of last year. Unimproved staggered fermions 
show taste breaking that is of the same order as the light pion masses, so that €?(a 2 ) ~ &{p ) is the 
appropriate power counting. Furthermore, the slopes differ significantly for the different tastes, in- 
dicating an important next-to-leading-order (NLO) contribution from &(a 2 p 2 ) terms. By contrast, 
for HYP-smeared fermions both the splittings and the differences in the slopes are much smaller. 
Thus we conclude that one should treat &(a 2 ) effects as being of NLO for HYP fermions on these 
lattices. The important practical question is whether this holds also for unquenched lattices. 

3. Comparison of asqtad and HYP-smeared staggered fermions 

To address this we have compared asqtad and HYP-smeared staggered valence fermions on 
one set of the coarse 2 + 1 flavor MILC lattices [Q]. The parameters of the study are given in 
Table ||. The range of quark masses for asqtad quarks is similar to that used above, the physical 
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Figure 3: (rim n ) 2 vs. 2r\m q for asqtad valence quarks on unquenched lattices. From Ref. 

strange quark mass am^ hys 0.045 lying near the upper end. It turns out, however, that the HYP 
quarks we use are somewhat lighter. For the same bare quark mass, the physical quark mass is 
about 1.4 times smaller with HYP than with asqtad quarks. Part of this difference is due to the 
normalization convention for tadpole improvement used by Ref. which has the effect that the 
bare quark mass for asqtad quarks should be multiplied by 1 / uq to match the convention used for 
unimproved and HYP quarks. Here uq is the "mean link", which is m 0.87 for these lattices. The 
remainder of the difference is presumably due to a difference in the values of Z m for the two lattices. 

In Fig. ||, we show the results obtained in Ref. [Q] using asqtad valence quarks. One sees 
that taste-breaking is significant, with ^(a 2 ) ~ &{p 2 ) being the appropriate power counting. Note 
that the taste-breaking with unimproved staggered quarks would be larger still, since these lattices 
are coarser than the quenched lattices used earlier. Despite the relatively large taste breaking, 



parameter 



value 



gauge action 
sea quarks 

a 

geometry 
# of confs 
valence quark type 
valence quark mass (asqtad) 
valence quark mass (HYP) 



1-loop tadpole-improved Symanzik 
2+1 flavors of asqtad staggered {ami = 0.01, am s = 0.05) 

6.76 
0.125 fm 
20 3 x 64 
640 (asqtad)/ 406 (HYP) 
asqtad and HYP staggered 
0.007, 0.01, 0.02, 0.03, 0.04, 0.05 
0.05, 0.01, 0.015, 0.02, 0.025, 0.03, 0.035, 0.04, 0.045, 0.05 



Table 2: Simulation parameters used for the comparison of asqtad and HYP valence quarks. Asqtad results 
are from Ref. [§. 
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Figure 4: {r\m % ) 2 vs. 2r\m q for HYP valence quarks on unquenched lattices, using cubic U(l) sources. Left 
panel: LT tastes; Right panel: NLT tastes. 



discretization effects of NLO (i.e. &{a 2 p 2 )) are very small. This is shown both by the lack of 
breaking of SO(4) symmetry, and the result that the slopes for different tastes are close. The latter 
result differs from that found above with unimproved fermions. 

Our new results using HYP-smeared staggered quarks are shown in Fig. ||. Results are with 
cubic U(l) sources; those with cubic wall sources are not yet available. We plot using the same 
dimensionless variables as in Fig. ||. The difference in slopes between asqtad and HYP fermions (by 
a factor of about 1.4) is due to the difference in quark mass conventions and renormalization factors 
discussed above. Irrespective of this detail, we see that taste breaking is substantially reduced for 
HYP-smeared quarks. To make this quantitative we quote, in Table [| the splittings in the chiral 
limit (using r\ to set the scale, rather than the 1/a used above). We see that HYP smearing reduces 
taste breaking by a factor of 2.5-3, although the ordering of states is the same for both asqtad and 
HYP actions. The improvement in taste-breaking is similar to that found using the HISQ action [g]. 



Taste (F) 


Af (asqtad) 


A F (HYP) 


Taste (F) 


Af (asqtad) 


Af (HYP) 




0.205(2) 


0.072(2) 




0.205(2) 


0.070(18) 




0.327(4) 


0.138(4) 




0.327(4) 


0.142(20) 




0.439(5) 


0.203(7) 


& 


0.439(5) 


0.210(28) 










0.537(15) 


0.259(39) 



Table 3: Taste symmetry breaking effect: Af = [r\m n (F)] 2 — [rim n (^s)] extrapolated to the chiral limit. 
Results from unquenched lattices comparing valence asqtad (from Ref. ^ — with SO(4) symmetry enforced) 
to HYP-smeared fermions (with cubic U(l) sources). 

We can characterize the taste-breaking errors in terms of a non-perturbative scale. Using 
the splitting of the "2-link" pions (taste ^ v ) from the Goldstone pion as an "average" splitting, 
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and setting Aa 2 m| = a 4 A 4 , we find A « 0.6 GeV. This should be compared to scales « 1 GeV 
for unimproved staggered fermions. Thus improvement has reduced discretization errors to the 
expected generic size, i.e. A ~ Aqcd- 

Looking in more detail at the results for HYP fermions, we note that SO(4) breaking is small, 
and that the slopes of m 2 n versus m graphs are similar for all tastes. Thus, as for the asqtad quarks, 
the NLO &(a 2 p 2 ) terms responsible for these effects must be small. 

Despite the reduction in taste breaking, it is clear from the figures that, for this range of quark 
masses (which is the range we are using for calculating matrix elements), the appropriate power 
counting for the MILC coarse lattices remains 0(a 2 ) ~ &{p 2 ) for HYP fermions. 

4. Conclusions 

The most important conclusion for future work is that HYP smearing leads to a significant 
reduction in taste breaking compared to the asqtad action. This was expected from previous studies, 
but it is important to demonstrate this on the unquenched MILC lattices that we are using for 
calculating weak matrix elements. The reduction is, however, not small enough to change the 
power counting for staggered chiral perturbation theory on the coarse lattices. If, however, one 
extrapolates to the fine MILC lattices, on which one expects that taste breaking will be reduced 
by about 2.5 (which is the relative size of a 2 a 2 ), then it seems likely that one will be able to treat 
the discretization errors as a NLO effect, i.e. that &{a 2 ) ~ 0{p ). Pictorially, the spectrum on the 
fine lattices should look like that for the quenched lattice results of Fig. ^[ This would lead to a 
significant simplification in the fitting forms from staggered chiral perturbation theory. 
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